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1.   INTRODUCTION 

Simulation  is  a  useful  tool  of  the  operations  analyst.  This  is 
not  a  new  concept,  the  first  recorded  simulations  were  conducted  by 
the  Chinese  in  the  form  of  "war  chess"  and  were  probably  used  to  teach 
young  men  some  of  the  concepts  of  battle  without  the  inherent  danger 
of  loss  of  life.  Later  accounts  of  simulations  have  been  recorded  by 
the  Prussians,  French  and  Germans  [7]. 

This  technique  of  analyzing  a  problem  by  simulation  is  now  employed 
by  all  branches  of  the  scientific  community.  The  present  day  high  speed 
digital  computer  has  given  rise  to  a  rapid  expansion  in  the  use  of 
simulation  as  a  method  of  solution  to  military,  scientific,  management 
and  many  other  types  of  problems. 

Convincing  the  reader  or  observer  that  a  simulation  "models  the 
real  world"  is  one  of  the  primary  problems  confronting  the  analyst 
who  uses  simulation  techniques  in  the  solution  of  a  problem.  One  way 
to  minimize  this  doubt  is  to  increase  the  role  of  the  human  in  the 
simulation.  This  can,  in  many  cases,  be  done  by  the  techniques  of  on- 
line simulation.  While  on  the  one  hand,  on-line  simulation  increases 
the  complexity  of  the  problem  by  nature  of  the  man  machine  interface, 
at  the  same  time,  on-line  simulation  adds  very  complex  logic  (the  man) 
to  the  problem  with  minimal  effort  on  the  part  of  the  designer. 

The  design  of  a  simulation,  and  in  particular  an  on-line  simulation, 
can  appear  to  be  a  formidable  task  when  first  considered.  However,  if 
the  designer  has  an  approach  in  mind  and  proceeds  in  an  organized 
manner,  the  problem  usually  divides  into  small  parts  that  are  each 
relatively  simple.   It  is  the  aim  of  this  thesis  to  take  a  representative 


problem  and  develop  a  computer  simulation  that  can  be  used  as  a 
reference  for  the  construction  of  on-line,  real  time,  computer  simu- 
lations in  general. 


2.  FORMULATION  OF  ON-LINE  SIMULATIONS 

In  the  formulation  of  on-line  simulations  the  designer  must  first 
lay  down  a  good  foundation  in  the  form  of  a  well-planned  outline.  This 
outline  must  begin  by  initializing  and  setting  the  scene  for  the  simu- 
lation. When  this  has  been  accomplished  the  designer  must  turn  his 
attention  to  the  development  of  a  loop  that  will  include  all  the  actions 
and  interactions  expected  to  occur  in  the  simulation.   Included  as  an 
integral  part  of  this  loop  is  a  timing  mechanism  that  is  flexible  enough 
to  allow  the  simulation  to  proceed  at  any  rate  required.  A  third  part 
of  the  formulation  is  concerned  with  providing  a  critique  of  the  simu- 
lation, either  as  a  running  critique  or  a  compilation  of  pertinent  facts 
at  the  end  of  the  simulation. 

INITIALIZING 

Initializing  is  the  term  applied  to  that  portion  of  the  simulation 
which  is  executed  before  "play"  begins.   It  includes  data  input,  assign- 
ment of  particular  values  to  the  parameters  and  entering  starting  values 
required  for  indexing  the  logic.  Herein  is  provided  the  flecibility 
required  of  any  useful  simulation.   In  the  initializing  portion  the  ground 
work  must  be  formed  for  performing  sensitivity  analysis  if  such  analysis 
is  required.  The  initializing  portion  of  the  simulation  must  allow  enough 
flexibility  to  provide  for  ti  e  various  scenarios  possible  in  the  particular 
simulation.  Clearly  tnen,  the  initializing  portion  must  be  designed  with 
these  purposes  as  primary  criteria  and  subject  to  boundary  conditions, 
such  as  equipment  capabilities. 


AN  ITERATIVE  LOOP 

In  general,  many  on-line  simulations  contain  an  iterative  loop  that 
is  cycled  for  each  time  period.  Therefore,  one  of  the  first  considerations 
to  be  made  in  this  phase  of  the  design  is  the  determination  of  the  stepping 
interval.   Several  factors  are  involved  in  this  selection,  the  most 
important  being  the  assurance  that  the  stepping  interval  is  compatible 
with  the  logic  flow  of  the  situation  being  simulated.  The  designer  must 
then  consider  the  amount  of  time  required  to  accomplish  the  most  com- 
plicated simulation  situation  that  can  occur  in  one  cycle.  These  con- 
siderations complete,  the  simulator  must  insure  that  the  cycling  is  such 
that  the  player  is  not  bored  with  the  data/action  as  presented  and  also 
that  this  data/action  is  not  presented  at  a  rate  too  rapid  for  the  player 
to  fully  comprehend.  With  these  considerations  in  mind,  a  tentative 
looping  cycle  can  be  constructed  and  the  designer  may  continue  to  develop 
the  necessary  routines  required  to  complete  the  iterative  loop. 

The  loop  should  now  be  fashioned  in  its  most  elementary  form.  The 
designer  must  consider  several  tasks  that  must  be  accomplished  during 
each  cycle.  These  include,  for  example,  advancing  all  participants  one 
time  cycle,  consideration  of  the  possible  interactions  that  can  occur  due 
to  these  moves,  tabulation  of  the  results  of  such  interactions,  presen- 
tation of  output  to  the  player,  permitting  the  player  to  communicate  with 
the  simulation,  delaying  the  next  cycle  until  the  proper  time  interval 
has  transpired,  and  possible  other  considerations  dependent  upon  the 
particular  simulation. 

Care  and  planning  must  be  exercised  in  the  construction  of  this  loop 
since  this  is  the  foundation  upon  which  the  designer  is  to  build  his 
simulation.   If  logical  errors  appear  in  the  order  of  these  routines  or 
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a  component  is  not  considered  in  this  loop,  the  remaining  portion  of 
the  design  will  be  difficult,  if  not  impossible.   Planning  in  this 
portion  of  the  development  will  be  time  well  spent. 

THE  CRITIQUE 

The  critique  is  that  portion  of  the  simulation  in  which  the  entire 
simulation,  or  any  integral  part,  is  analyzed  and  the  results  compiled 
in  condensed  form.  There  are  several  basic  tecnniques  that  may  be 
explored  in  this  part  of  the  simulation.  The  least  complex  of  which  is, 
in  most  cases,  a  complete  "recording"  of  the  game  that  can  be  "replayed" 
at  a  later  time  at  any  speed  desired.  A  more  complex  approach  to  the 
problem  of  performing  a  critique  of  a  simulation  is  that  of  including 
a  recording/analysing  routine  in  the  iterative  loop.  This  routine  would 
extract  the  desired  information  during  each  cycle  of  the  loop.   It  is 
apparent  that  not  every  cycle  need  be  recorded,  and  therefore,  a  decision 
logic  must  be  included  that  will  extract  all  necessary  information. 
Associated  with  this  technique  must  be  a  recording  routine  that  can  be 
queried  at  intervals  or  at  the  conclusion  of  the  simulation. 

Much  care  must  be  exercised  in  formulating  the  design  and  location 
of  the  critique  routine  since  this  is  the  major  mode  of  analysis  available 
to  the  analyst.  Flexibility  and  adaptability  are  considerations  that 
must  be  made  to  allow  che  simulation  to  be  fully  appreciated  as  an 
analytical  tool  and  not  just  a  "parlor  game." 
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3.   PROGRAM  MODULARIZATION 

In  general,  programmers  attempt  to  modularize  their  programs.  There 
are  several  reasons  for  this,  the  most  obvious  being  that  of  providing 
logical  grouping  of  ideas.  The  large  scale  simulation  is  usually  modu- 
larized by  the  use  of  subroutines.  This  use  of  subroutines  is  convenient 
because: 

a.  Several  separate  groups  may  be  working  on  various  sections  of 
the  problem,  and  in  many  instances  the  use  of  subroutines  is  the  best 
technique. 

b.  Computations  which  are  to  be  called  upon  several  times  in  the 
main  program  are  best  handled  by  the  use  of  subroutines. 

c.  The  program  may  be  of  a  magnitude  such  that  the  entire  program 
cannot  be  compiled  in  one  pass. 

These  reasons  are  valid  for  the  large  scale  simulation.   In  simu- 
lations which  are  moderate  to  small  in  size  the  use  of  subroutines  may 
add  unnecessary  factors  to  be  considered,  with  the  exception  of  reason 
(b)  above  which  is  a  valid  reason  for  the  use  of  subroutines  in  most 
computer  applications. 

By  careful  construction  of  the  statement  numbering  scheme,  available 
in  languages  such  as  FORTRAN,  the  programmer  of  the  moderate  to  small  size 
simulation  can  modularize  his  program  without  the  additional  consideration 
of  designating  common  storage  and  the  other  difficulties  experienced  when 
programming  subroutines. 

ADVANTAGES 
The  above  technique  allows  the  designer  to  use  the  familiar  computer 
languages,  such  as  FORTRAN,  in  place  of  special  simulation  languages. 
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It  also  alleviates  the  requirement  of  either  providing  the  reader  with  a 
description  of  a  special  language  or  causing  him  to  go  to  another  source 
to  interpret  the  program.   Using  a  "standard"  language,  such  as  FORTRAN, 
the  designer  can  reasonably  assume  that  the  reader  needs  little  or  no 
explanation. 

FLEXIBILITY 
The  use  of  modularization  with  statement  numbers  gives  all  the 
inherent  flexibility  observed  in  the  special  languages  when  the  simulation 
is  of  such  a  magnitude  as  to  allow  compilation  in  one  pass.  When  simu- 
lations are  in  the  design  and  programming  stages,  the  use  of  the  above 
technique  can  allow  as  much  flexibility  as  the  simulation  requires.   Exit 
from  the  blocks  when  programming  in  FORTRAN  can  be  accomplished  at  any 
logical  point  with  nothing  more  than  a  simple  GO  TO  or  COMPUTED  GO  TO 
statement,  and  thus  the  logic  flow  is  easily  accomplished  by  this  method. 
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4.   SYSTEM  VARIABLES,  COORDINATES  AND  NOTATION 

Certain  special  simulation  languages  provide  for  dynamic  allocation 
of  storage  for  tables  which  facilitates  the  designation  of  data  in  a 
flexible  and  expandable  form.  This  is  not  necessarily  required  by  the 
moderate  size  simulation.  The  designer  of  these  moderate  size  simulations 
may  find  the  techniques  as  explained  below  more  desirable,  since  his 
problem  is  not  generally  one  of  storage  limitations.  The  criterion  for 
the  dynamic  versus  preset  storage  decision  is  felt  to  be  that  of  program 
size.   In  programming  the  moderate  size  simulation  the  programmer  may 
find  that  constructing  a  scheme  for  naming  variables  may  create  easier  to 
work  with  variables  than  the  complex  tabular  form  of  the  special  simu- 
lation language.  This  was  found  to  be  true  in  the  programming  of  the 
example  simulation,  MULNUCl,  of  this  thesis. 

The  general  approach  to  the  naming  of  variables  in  the  example 
program  was  that  of  using  vectors  to  represent  a  given  parameter,  each 
component  of  the  vector  being  representative  of  the  value  of  that  parameter 
with  respect  to  the  unit  concerned.  An  example,  that  of  the  X-coordinate 
of  the  i*""  destroyer,  being  DDX(I)  .   Once  the  scheme  is  understood  the 
programming  moves  along  without  a  great  amount  of  thought  required  by  the 
programmer  as  far  as  variable  names  are  concerned.  Certain  variables 
must  be  set  aside  as  dummy  or  temporary  and  these  logically  take  on  forms 
such  as:  ITEMP,  TEMP3,  DUMMY ( I) ,  and  forms  that  immediately  classify  them 
in  this  category. 

In  general,  several  different  coordinate  systems  are  required  in 
simulations.   In  the  war  game  an  overall  "area  of  play"  must  be  established, 
This  can  be  either  rectangular  or  polar,  two  or  three  dimensional.  The 
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axes  of  the  rectangular  coordinate  system  are  not  necessarily  graduated 
in  the  same  units.   In  the  example  program,  for  instance,  it  will  be 
seen  that  the  third  axis  (depth)  is  dimensioned  in  feet,  while  the  two 
major  axes  are  dimensioned  in  yards.  Further,  the  third  axis  has  the 
normally  negative  direction  established  as  positive.  These,  perhaps 
unorthodox,  measures  are  taken  to  facilitate  programming;  however,  they 
must  be  spelled  out  in  the  documentation  of  the  simulation  to  prevent 
possible  misunderstanding. 

Many  simulations  require  more  than  one  coordinate  system  to  be 
employed.  The  overall  play  is  perhaps  in  a  rectangular  coordinate  system, 
while  range  and  bearing  information  may  be  required  during  the  play  of 
the  game.  This  will  necessitate  the  incorporation  of  an  overlay  of  one, 
or  perhaps  several,  polar  coordinate  systems  upon  the  base  system.   If 
a  "close-up"  view  is  required  during  the  play  a  translation  and/or 
expansion  to  another  rectangular  system  may  be  required.   It  can  now  be 
seen  that,  in  general,  several  cootdinate  systems  will  be  used  in  a  war 
game  type  of  simulation.  Therefore,  a  plan  for  the  designation  of  these 
various  coordinate  systems  and  their  respective  transformations  must  be 
established  early  in  the  formulation  of  a  simulation. 
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5.   MULNUC1  -  AN  EXAMPLE 

MULNUCl  is  an  on-line,  real  time  simulation  used  as  an  example  of 
the  application  of  ideas  presented  in  sections  one  through  four.  The 
simulation  had  its  beginning  at  the  Naval  Radiological  Defense  Laboratory, 
Hunters  Point,  in  the  summer  of  1965.  During  a  six-week  tour  of  the 
Naval  Radiological  Defense  Laboratory,  it  was  found  that  little  had  been 
done  in  the  exploration  of  tactics  and  possible  reactions  of  surface 
antisubmarine  destroyers  exposed  to  a  self-inflicted  multiple  burst 
nuclear  environment.  At  this  time  the  simulation  used  as  an  example  in 
this  thesis  had  its  beginning.   In  the  example  program,  MULNUCl,  all 
classified  input  parameters  have  been  assigned  fictitious  values  so  that 
the  computer  program,  as  presented  in  this  thesis,  could  remain  unclassified. 

In  1965,  Lieutenant  J.  E.  Johnson  programmed  the  on-line  display, 
Display  Data  Corporation  model  DD  65,  using  a  rather  simple  simulation 
situation  [2J.  He  did  make  a  contribution  in  the  form  of  an  advancement 
in  the  techniques  of  on-line  programming  of  simulations.  Having  observed 
a  demonstration  of  Johnson's  program,  it  was  felt  that  the  technique  of 
on-line  display  would  be  ideal  for  the  envisioned  program,  MULNUCl. 

Several  links  were  missing  in  the  chain  necessary  to  put  the  envisioned 
program  on-line.  The  first  link  was  a  requirement  for  a  routine  to  gen- 
erate circles  of  arbitrary  size  and  location.  This  was  accomplished  with 
subroutine  CIRCLE  (see  Appendix  III).  After  completing  subroutine 
CIRCLE,  attention  was  turned  to  the  necessity  for  a  random  number  gen- 
erator capable  of  generating  several  types  of  random  variables.  The 
distributions  of  random  variables  required  were  uniform,  normal,  and 
circular  normal.  These  generators  were  written  in  the  form  of  subroutines 
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UNIFORM,  NORMAL,  and  ERROR.  The  subroutine  RANVAR  is  the  basic  random 
number  generator  called  by  these  subroutines  in  the  generation  of  their 
respective  random  variables.  The  above  routines  completed,  only  the 
communication  routines  required  to  link  the  Control  Data  Corporation 
1604  and  160  computers  remained.  This  requirement  was  satisfied  by  sub- 
routines DCIRCLE,  DTRACK,  PARAMS  and  DSTATUS  (see  Appendix  III  and 
Acknowledgements) . 

At  this  point  the  preliminary  work  was  complete  and  the  formulation 
of  the  initializing,  iterative  loop,  and  critique  portions  of  the  simu- 
lation was  begun.  These  three  basic  steps,  as  discussed  in  section  two, 
were  incorporated  into  an  executive  control  block. 

EXECUTIVE  CONTROL 
The  Executive  Control  routine  is  flow  charted  in  Appendix  II  and 
consists  of  three  major  parts.  The  first  of  these  parts  is  the 
initializing  portion.   It  is  made  up  of  four  blocks: 

1.  Inputs 

2.  Set  Constants 

3.  Initialize 

4.  Enter  Input  Changes 

In  this  subsection  we  shall  consider  the  first  three,  leaving  the 
latter  for  discussion  in  the  subsection  titled  Man  Machine  Interface. 
The  Inputs  block  is  the  one  in  which  "standard"  or  nominal  input 
parameters  are  set.  These  parameters  are  listed  below. 
Number  of  Destroyers 
Time  Factor 

Destroyer  and  Submarine  Maximum  Speed 
Submarine  Hull  Parameters 
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Initial  Positions 

Initial  Courses  and  Speeds 

Nominal  Yield  of  Nuclear  Weapons 

Depth  of  Thermocline 

Maximum  Range  of  Weapons 

Detonation  Parameters  for  Weapons 

Random  Number  Generator  Initializer 
Any  of  these  parameters  may  be  changed  in  the  enter  input  changes  block 
(see  Man  Machine  Interface  subsection) .   These  parameters  were  chosen 
as  the  minimal  requirements  necessary  to  produce  a  simulation  that  has 
some  realism  and  yet  is  not  too  complex.  The  structure  of  this  program 
is  such  that  any  block  can  be  expanded  to  include  more  parameters , 
thereby  creating  a  more  realistic  simulation. 

Set  Constants  in  a  block  used,  as  the  name  implies,  to  initialize 
non-changeable  inputs.   In  this  block,  all  the  logic  indicators  are  set 
to  orient  the  game.  All  damage  and  radiation  levels  are  set  at  zero. 
The  indices  for  tracking,  firing,  and  sonar  contact  are  set  at  zero. 
This  is  easily  followed  by  cross  referencing  Appendices  I  and  IV. 

The  Initializing  block  begins  by  initializing  the  random  number 
generator  subroutine  and  then  calculates  the  following  parameters: 

Water  Temperature  Gradient 

Submarine  Crush  Depth 

Operational  Depth  of  Submarine 

Wind  Direction  and  Velocity 

Minimum  Safe  Range  of  Weapons  [6] 

Effective  Sonar  Range 

Sea  State 
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These  calculations  are  straight  forward  and  can  easily  be  followed  by 
cross  referencing  Appendices  I  and  IV.  The  block  also  presents  input 
parameters  of  interest  to  the  destroyer  team,  submarine  team  (if  selected) , 
and  critique  routine. 

Subsurface  nuclear  bursts  are  divided  into  four  classifications: 

1.  Very  Shallow 

2.  Shallow 

3.  Deep 

4.  Very  Deep  [4] 

The  criteria  for  selection  of  classification  are  depth  of  burst  and  yield. 
The  determination  of  classification  of  burst  is  made  at  this  point  in  the 
program  and  an  index  IDEEP  is  set  (see  Appendix  I) .  The  four  matrices  of 
output  data  are  filled  with  negative  zero,  since  negative  zero  is  pro- 
grammed not  to  print  on  the  display.  Finally,  if  the  role  of  the  sub- 
marine is  to  be  played  by  the  computer,  the  basic  strategy  of  the  sub- 
marine is  randomly  determined  (see  Submarine  Logic  Model) . 

The  next  major  part  of  the  Executive  Control  Routine  is  the  iterative 
loop.   It  is  made  up  of  eight  blocks: 

1.  Plot  Positions 

2.  Display  Data 

3.  Plot  Generator 

4.  Interactions 

5.  Radiation  Model 

6.  Enter  Changes 

7.  Submarine  Logic  Model 

8.  Time  Loop 

The  first  three  of  these  will  be  considered  now,  while  the  remaining  blocks 
will  be  explained  in  later  subsections. 
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The  Plot  Positions  block  is  utilized  to  transmit  output  information 
to  the  left  tube  of  the  DD  65  (see  Figure  2).  This  output  information 
consists  of  the  following  data: 

Destroyer  Tracks 

Sonar  Contact  Plots 

Destroyer  Courses,  Speeds  and  Coordinates 

Orientation  and  Size  of  Area  Displayed 

Important  Messages  to  the  Player 
This  information  presentation  is  covered  in  more  detail  in  the  section 
on  Man  Machine  Interface* 

The  Display  Data  block  performs  the  same  function  with  respect  to 
the  right  tube  of  the  DD  65  (see  Figure  3).  Figure  3  lists  the  data 
displayed  by  this  block  and  for  this  reason  it  will  not  be  listed  at 
this  point. 

The  final  major  part  of  the  Executive  Control  Routine  to  be  con- 
sidered at  this  time  is  the  Critique  block.  This  block  critiques  the 
simulation  by  performing  several  tasks.  The  entire  simulation  is  recorded 
on  magnetic  tape  and  can  be  reviewed  at  a  later  time.  If  any  significant 
action  or  interaction  occurs  game  time  and  the  nature  of  the  action/ 
interaction  are  recorded  by  the  Critique  I  block.  A  narrative  print  out 
is  made  at  the  conclusion  of  the  simulation  by  the  Critique  II  block. 
To  correlate  this  information  a  graph  plot  of  the  tracks  of  destroyers, 
submarines,  and  all  pools  and  clouds  of  radiation  is  made.  Detail  of 
these  operations  is  documented  in  Appendix  IV. 

PLOT  GENERATOR 

The  Plot  Generator  block  advances  all  participants  each  time  step 

(see  Appendix  II).  The  block  then  determines  if  there  are  any  clouds  or 

pools  of  radiation.   If  there  are  clouds,  they  are  advanced  using  wind 
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Rt^HT    TUBE     DtSPlAX 


INPUT  PARAMEIERS 

NUMBER  OF  DESTROYERS  ...   2 

SIZE  OF  ASROC  WARHEAD  ...  2.0   KT 

DEPTH  OF  THENMOCLINE  ...120.0   FT 

MAXIMUM  RANGE  OF  ASROC  .  9000.0  YARDS 

DEPTH  OF  BURST 700.0  FT 

SINK  RATE  OF  WARHEAD  ...  18.0    FT/SEC 

TIME  FACTOR  5.0 

WATER    TEMP    GRADIENT    ....    -.335    DEG/100FT 


^FF    SONAR  RANGE 

6171. 

*IND  DIRECTION 

179 

GAME  TIME 

42.5 

WIND  VELOCITY 

28. 

WARHEAD  SIZE 

2.0 

SEA  STATE 

7 

DD1 

DD2 

MAX  SPEED  AVAIL 

33. 

33. 

CONTACT  CLASS 

1 

SONAR  RANGE 

4899. 

SONAR  BEARING 

65 

DOPPLER 

1 

TARGET  COURSE 

2 

TARGET  SPEED 

22. 

FIRING  SOLU1 ION 

1 

RADIATION  RATE 

• 

• 

RADIATION  DOSE 

• 

. 

FIGURE   3 


velocity  and  direction  to  determine  motion,  and  the  radius  is  computed 

[1],   Pools  are  considered  to  be  stationary  since  their  drift  is  negligible. 

INTERACTIONS 
Interaction  is  the  largest  block  and  for  clarity  has  been  broken 
into  sub-blocks.  These  sub-blocks  are: 

1.  Sonar  Contact  Model 

2.  Contact  Tracking  Model 

3.  Weapon  Firing  Model 

4.  Evaluation  Model 

The  interactions  block  is  constructed  in  the  form  of  a  loop  that  considers 
the  interactions  of  each  unit  in  succession.  The  organization  is  well 
documented  in  Appendices  II  and  IV  and  will  not  be  further  explored  at 
this  time,  however,  the  details  of  the  four  sub-blocks  listed  will  be 
considered  below. 

SONAR  CONTACT  MODEL 
The  Sonar  Contact  Model  uses  a  ray  path  theory  detection  scheme  in 
a  deterministic  manner  [3].  This  deterministic  detection  range  then  has 
a  variance  superimposed  upon  it.  The  net  result  is  a  fairly  realistic 
sonar  detection  model.  The  major  limitation  of  this  model  is  that  it 
only  handles  the  constant  temperature  gradient  case.  The  inclusion  of 
other  gradients  causes  the  sonar  detection  problem  to  assume  a  much  more 
complex  nature.   Included  in  the  model  is  a  degradation  of  sonar  range 
due  to  excessive  destroyer  speed.  The  range  and  bearing  given  as  outputs 
from  this  model  have  range  and  bearing  errors  included.  The  data,  with 
these  errors,  is  then  utilized  by  the  tracking  and  firing  models.  Doppler 
is  also  calculated  in  the  model  and  sent  to  the  display  as: 
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1.  Up  Doppler 

2.  Down  Doppler 
3*  No  Doppler 

CONTACT  TRACKING  MODEL 
The  Contact  Tracking  Model  uses  as  input  data  the  output  of  the 
Sonar  Contact  Model.  A  simple  criteria,  requiring  three  consecutive 
marks  from  the  sonar  model,  is  used  to  distinguish  non-contacts  from 
contacts.  Once  three  consecutive  marks  are  received,  the  model  determines 
the  contact  course  and  speed.  This  is  done  with  a  simple  no  parameter 
track  model  that  considers  at  least  three  but  no  more  than  five  marks. 
The  course  and  speed  of  the  contact  are  determined  from  the  first  and 
last  of  these  marks.  This  track  model  is  unsophisticated  but  the  error 
induced  in  the  output  is  fairly  realistic.  This  model  is  of  the  first 
that  should  be  improved  upon  if  more  work  is  to  be  done  on  this  simulation. 
The  output  of  this  model  is  in  the  form  of  contact  course  and  speed.  This 
model  provides  dual  routing,  dependent  upon  the  track.  This  will  be 
discussed  further  in  the  subsection  on  Block  Sequencing.  Determination 
of  whether  or  not  the  target  is  in  firing  range  is  made  just  before  exiting 
the  routine. 

WEAPON  FIRING  MODEL 
The  Weapon  Firing  Model  takes  the  last  position  of  the  contact  from 
the  Sonar  Model,  the  contact  course  and  speed  from  the  Tracking  Model, 
and  then  determines  a  firing  solution.  The  range  is  considered,  with 
time  of  flight  and  sink  time,  and  the  time  of  burst  is  determined.  A 
dead  reckoning  position  of  the  target  is  computed  from  the  track  data 
and  this  position  becomes  the  aim  point  of  the  weapon.  The  model  then 
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calls  subroutine  ERROR  from  which  the  true  fall  of  shot  is  determined. 
The  weapon  is  given  a  reliability  check  in  the  model  and  if  this  test  is 
failed,  the  player  will  be  notified  that  the  weapon  has  misfired  (see 
Man  Machine  Interface  subsection) . 

EVALUATION  MODEL 

The  evaluation  model  initializes  the  pool  and  cloud  of  radiation 
created  by  the  subsurface  nuclear  burst.  To  accomplish  this  the  model 
takes  data  from  the  firing  model  for  the  location  of  ground  zero  and 
data  from  the  inputs  block  for  yield,  depth  of  burst,  and  type  of  burst. 
The  radius  of  the  radioactive  pool  is  then  determined  [4]  (the  cloud 
parameters  are  computed  in  the  plot  generator  model) . 

A  simple  criteria  for  damage  to  the  submarine  is  used.  The  lethal 
range  is  determined  using  submarine  hull  parameters,  yield  of  warhead, 
and  submarine  depth  as  received  from  the  submarine  logic  model.   Slant 
range  to  the  burst  from  the  submarine  is  computed.  If  the  submarine  is 
within  the  lethal  range,  damage  is  100%.  If  the  submarine  is  outside  a 
radius  equal  to  twice  the  lethal  range,  damage  is  zero.  Values  of  sub- 
marine damage  between  zero  and  100%  are  computed  by  a  linear  relation- 
ship, then,  if  at  any  time  the  submarine's  total  damage  reaches  the  75% 
level,  the  game  is  terminated  with  the  submarine  considered  as  having 
been  sunk. 

RADIATION  MODEL 
The  Radiation  Model  determines  if  any  weapons  have  been  detonated. 
If  none  have  been  detonated,  the  block  is  bypassed  and  the  program  con- 
tinues.  If  weapons  have  been  detonated,  the  model  computes  the  distance 
of  each  destroyer  from  all  pools  and  clouds  of  radiation.  The  location 
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and  size  of  all  pools  and  clouds  is  received  from  the  plot  generator 
model  along  with  the  location  of  the  destroyers.  The  model  then  determines 
if  the  destroyers  are  within  the  perimeter  of  any  pool  or  cloud.  If  thi6 
condition  exists,  the  radiation  level  in  the  cloud  or  pool  is  calculated 
[1,  4].  The  total  radiation  being  received  by  each  destroyer  is  then 
calculated.  The  radiation  rate  and  total  radiation  dose  for  each 
destroyer  are  computed  and  sent  to  the  display  (see  Man  Machine  Interface 
subsection)  . 

SUBMARINE  LOGIC  MODEL 
The  submarine  can  be  controlled  in  two  ways: 

1.  By  a  submarine  team. 

2.  By  the  computer. 

This  decision  is  made  in  the  enter  input  changes  routine.  The  program 
is  such  that  the  computer  will  play  the  role  of  the  submarine  unless  the 
variable  ISUB  is  set  equal  to  one  by  the  enter  input  changes  routine. 
If  ISUB  is  set  equal  to  one,  control  of  submarine  depth,  course,  speed, 
and  weapon  firing  is  turned  over  to  the  submarine  team.  This  team  will 
receive  passive  sonar  bearings  and  screw  beat  information  from  the 
console  typewriter  of  the  CDC  1604.   They  will  be  able  to  control  the 
movements  of  the  submarine  and  its  weapon  firing  by  means  of  the  console 
typewriter  and  selective  jump  keys.  The  program  will  query  selective 
jump  key  number  two,  once  each  cycle,  to  determine  if  orders  to  the 
submarine  are  to  be  received.  If  selective  jump  key  two  is  set  the 
computer  will  request  orders  (see  Man  Machine  Interface  subsection) . 
The  weapons  (torpedoes)  can  be  fired  by  setting  selective  jump  key  three 
on  the  CDC  1604  console  (see  Man  Machine  Interface  subsection) .  Under 
this  condition  the  simulation  becomes  a  conflict  between  two  teams: 
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1.  A  submarine  team  using  the  CDC  1604  console. 

2.  A  destroyer  team  using  the  DD  65  display  console. 

If  ISUB  is  unchanged  by  enter  input  changes  (inputs  block  sets 
this  variable  equal  to  zero) ,  the  moves  of  the  submarine  are  controlled 
by  the  computer.  This  being  the  case,  two  basic  initial  tactical  sit- 
uations are  available  to  the  player.  The  first  of  these  places  the 
submarine  on  the  surface,  at  the  origin  of  the  playing  area.  The  sub- 
marine knows  that  he  has  been  sighted  and  the  game  proceeds.  The  second 
situation  has  the  submarine  randomly  located  in  the  upper  half  of  the 
playing  area.   In  this  case  the  submarine's  position  is  not  known  by  the 
player.   It  should  be  noted  that  this  is  the  initial  situation  if  t?ie 
submarine  is  to  be  controlled  by  a  submarine  team. 

Initial  situation  one  is  selected  by  setting  INITIAL  equal  to  zero, 
while  situation  two  is  selected  by  setting  this  variable  equal  to  one. 
Having  selected  the  initial  situation  the  computer  then  selects  one 
of  three  basic  strategies: 

1.  The  submarine  runs  for  it. 

2.  The  submarine  tries  to  transit  between  the  two  destroyers. 

3.  The  submarine  tries  an  end  run,  flanking  the  two  destroyers. 
These  strategies  are  easily  followed  in  Appendix  II  and  will  not  be 
explored  further  at  this  point. 

SEQUENCING  OF  PROGRAM  BLOCKS 
The  blocks  in  this  simulation  are  of  two  basic  types: 

1.  A  single  point  of  exit • 

2.  Multiple  points  of  exit. 

The  block  that  has  a  single  exit  point  might  be  called  a  standard  block, 
in  that,  the  block  is  called  upon  to  perform  a  computation  but  no 
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branching  of  logic  is  done  within  the  block.  The  multiple  exit  block 
is  one  in  which  the  program  is  routed  differently  depending  upon  logic 
decisions  made  within  the  block.  The  sonar  contact  block  is  an  example 
of  this  type.   In  this  block  the  routing  depends  upon  the  results  of  a 
sonar  search.   If  contact  is  made,  the  block  exists  to  the  contact 
tracking  block.   If  no  contact  is  made,  the  block  exits  to  consider  the 
next  destroyer. 

The  executive  control  routine  flow  chart  in  Appendix  II  illustrates 
the  time  sequencing  of  the  major  program  blocks.  The  simulation  is 
delayed  at  four  points  in  the  program  which  are: 

1.  Enter  input  changes  routine. 

2.  Enter  changes. 

3.  Submarine  logic  block. 

4.  The  time  loop. 

The  first  of  these  interruptions  takes  place  only  during  the  initializing 
portion.  At  this  point  any  change  to  the  in^ut  parameters  is  made.  The 
second  of  these  interruptions  is  made  once  during  each  time  step.  This 
is  the  point  at  which  destroyer  team  changes  are  sent  to  the  CDC  1604. 
The  third  of  these  interruptions  takes  place  if  the  situation  using  the 
submarine  team  has  been  selected.   In  this  case  the  simulation  may  be 
interrupted  every  time  step  to  allow  the  submarine  team  to  enter  changes. 
The  fourth  of  these  interruptions  occurs  each  cycle  and  maintains  the 
time  stepping  interval. 

It  will  be  noted  that  in  the  executive  control  routine,  each  block 
is  considered  in  turn,  no  block  is  bypassed.  In  the  interactions  block, 
it  will  be  noted  that,  sub-blocks  are  not  always  considered.   No  sonar 
contact  by  the  sonar  contact  block  causes  the  contact  tracking  block  to 
be  bypassed.  The  same  is  true  of  the  weapon  firing  block  if  the  tracking 
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block  does  not  generate  a  satisfactory  track. 

MAN  MACHINE  INTERFACE 
This  subsection  is  concerned  with  communications,  both  into  and 
out  of  the  computer.  These  man-machine  communications  fall  into  four 
types: 

1.  Data  to  and  from  the  destroyer  team. 

2.  Data  to  and  from  the  submarine  team. 

3.  Data  to  the  various  modes  of  the  critique  routine. 

4.  Input  changes. 

The  first  type  breaks  into  three  parts: 

1.  Right  tube  information. 

2.  Left  tube  information. 

3.  Changes  sent  to  the  CDC  1604. 

The  right  tube  gives  the  destroyer  team  data  in  tabular  form  as  illustrated 
in  Figure  3.  The  left  tube  will  display  the  tracks  of  the  destroyers, 
any  sonar  contacts,  and  all  pools  and  clouds  of  radiation.  Also  included 
in  the  display  on  the  left  tube  is  a  series  of  windows  in  which  data  can 
be  displayed.  The  windows  are  numbered  as  shown  in  Figure  2.  Window 
data  assignments  are  as  listed  below. 

Window  Data  Assignment  Units 

1  X-coordinate  of  destroyer  1  100  yards 

2  X-coordinate  of  destroyer  2  100  yards 

3  Y-coordinate  of  destroyer  1  100  yards 

4  Y-coordinate  of  destroyer  2  100  yards 

5  Course  of  destroyer  1  Degrees  true 

6  Speed  of  destroyer  1  Knots 
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Window  Date  Assignment  Units 

7  Course  of  destroyer  2  Degrees  true 

8  Speed  of  destroyer  2  Knots 

9  X-coordinate  of  left  tube  center  100  yards 
10  Y-coordinate  of  left  tube  center  100  yards 
11-14  Available  for  flash  messages             Alfa-numeric 

15  Not  used 

16  Radius  of  display  on  left  tube  1000  yards 

Windows  5-10  and  16  are  controllable  from  the  DD  65  console  by  means  of 
a  discrete  digital  type  control  system.  In  this  manner  the  player  is  able 
to  change  the  destroyer  course  and  speed,  or  "zoom"  the  display  in  on  any 
point  in  the  area  of  play.  Let  us  first  consider  windows  five  through 
eight. 

On  the  DD  65  console  (see  Figure  4)  there  are  buttons  labelled 
CONN  DD  1  and  CONN  DD  2.  By  depressing  one  or  both  of  these  buttons  the 
player  is  given  control  of  the  course  and  speed  of  the  destroyer  or 
destroyers  selected.  Near  the  button  just  selected  is  a  group  of  four 
buttons  (see  Figure  4)  labelled  RIGHT,  LEFT,  FAST/UP,  and  SLOW/DOWN. 
Depressing  one  of  these  buttons  will  cause  the  appropriate  variable  in 
windows  five  through  eight  to  change  in  the  desired  direction.  As  an 
example,  if  the  player  depresses  both  CONN  buttons  and  then  the  button 
labelled  RIGHT  -  both  destroyers  will  commence  a  turn  to  the  right. 
They  will  continue  this  turn  until  the  buttcn  labelled  RIGHT  is  released. 
The  player  can  "come  to  a  course"  by  depressing  the  correct  button  until 
the  desired  course  is  displayed  in  the  respective  window.  The  same 
procedure  is  used  for  changing  speed.  It  should  be  noted  that  the 
windows  can  be  changed  to  values  that  are  unacceptable,  such  as  -5  knots, 
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in  this  case  the  program  will  cause  this  value  to  be  changed  back  into 
the  acceptable  range  during  the  following  cycle.  The  value  of  course 
can  be  increased  to  values  greater  than  360  degrees,  in  which  case  the 
program  will  correct  to  the  acceptable  valus.  For  example,  if  the  player 
increases  the  course  to  390  degrees,  the  program  will  convert  this  to 
030  degrees  during  the  next  cycle. 

The  same  basic  procedure  is  utilized  to  "zoom"  the  display  to  any 
location  desired.  First  the  player  selects  the  button  labelled  SHIFT 
(see  Figure  4).  He  now  has  control  of  windows  nine  and  ten.  Depressing 
the  button  labelled  RIGHT  will  cause  the  display  to  shift  to  the  right, 
LEFT  accomplishes  the  same  action  but  to  the  left,  and  similarly  with  UP 
and  DOWN.  To  zoom  in,  the  rotary  switch  in  the  upper  right  hand  corner 
of  the  console  is  used  (see  Figure  4).  This  switch  has  six  positions 
labelled  4,  8,  16,  32,  64,  and  128.  Selection  of  one  of  the  six  positions 
will  cause  the  radius  of  the  displayed  area  to  be  that  of  the  value 
selected,  in  thousands  of  yards.  As  an  example,  selecting  SHIFT  and 
changing  windows  nine  and  ten  to  240  and  150  respectively,  then  setting 
the  rotary  switch  to  eight,  will  cause  the  area  centered  at  (24000,15000) 
with  a  radius  of  8000  yards  to  be  displayed  on  the  left  tube.  Note  that 
window  16  will  show  the  value  eight,  while  windows  nine  and  ten  will  show 
240  and  150  respectively. 

Windows  11  through  14  are  used  to  send  the  player  the  following 

flash  messages. 

SONAR  CONTACT 

ASROC  FIRED 

ASROC  MISFIRE 

SUB  SUNK 

DD  SUNK 
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TARGET  IN  RANGE 
TARGET  TOO  CLOSE  TO  SHOOT 
TARGET  OUT  OF  RANGE 
Another  method  of  communications  is  available  to  the  destroyer  team. 
By  depressing  the  DD  1  FIRE  button,  destroyer  number  one  fires  an  ASROC. 
The  same  procedure  is  used  with  the  DD  2  FIRE  button. 

The  next  type  of  communications  available  is  that  of  the  submarine 
team.  This  is,  of  course,  non-existent  if  the  option  is  chosen  in  which 
the  submarine  is  played  by  the  computer.  The  submarine  team  will  receive 
messages  each  cycle  containing  bearing  and  screw  beat  information.  The 
submarine  team  may  then  choose  to  maneuver  the  submarine  by  setting 
selective  jump  key  number  two  on  the  console  of  CDC  1604.  The  computer 
will  then  type  COURSE  ORDERS  on  the  console  typewriter.  This  is  the 
indication  that  the  computer  is  ready  to  receive  course  changes.   If  no 
change  is  desired,  the  old  course  is  typed  in.   If  a  change  is  requested, 
the  new  course  is  typed  in.  The  course  typed  in  should  be  of  the  form 
090.  followed  by  a  carriage  return.  This  will  change  the  course  and  the 
computer  will  return  NEW  SUB  COURSE  090.  This  completes  the  course  change 
cycle  and  the  computer  will  then  type  SPEED  ORDERS,  the  same  procedure 
is  used  to  enter  speed  changes.  Upon  completion  of  the  speed  entry  the 
computer  will  return  NEW  SUB  SPEED  15,  followed  by  DEPTH  ORDERS.  The  new 
depth  is  now  entered,  and  the  computer  will  return  NEW  SUB  DEPTH  1050. 
The  routine  is  now  finished  and  the  program  continues. 

The  only  other  action  available  to  the  submarine  team  is  that  of 
firing  torpedoes.  This  is  accomplished  by  depressing  selective  jump  key 
number  three  until  the  typewriter  returns  TORPEDO  FIRED.  At  this  time 
the  selective  jump  key  should  be  returned  to  the  normal  position,  unless 
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another  torpedo  is  desired.  If  a  hit  is  scored,  the  typewriter  will 
return  DESTROYER  SUNK.  When  both  destroyers  are  sunk  the  typewriter 
will  return  GAME  OVER.  At  the  beginning  of  the  game,  the  typewriter 
will  give  the  submarine  team  the  following  information: 

1.  Maximum  speed  available  to  submarine. 

2.  Maximum  depth  allowable. 

The  third  type  of  communication  is  with  the  various  critique  routines. 
Critique  is  accomplished  in  three  ways: 

1.  A  recording  of  all  information  on  the  DD  65  display  is  recorded 
by  the  tape  unit  near  the  DD  65. 

2.  A  graph  of  the  tracks  of  the  DD's,  submarine,  and  all  radiation 
is  made  on  tape  unit  eight  of  the  CDC  1604  (the  graph  has  game  time 
recorded  by  each  mark  to  aid  in  correlating  with  the  various  other  parts 
of  the  critique  routine). 

3.  A  critique  of  all  important  events  and  their  time  is  recorded 

on  tape  unit  five  of  the  CDC  1604  for  print  out  at  the  conclusions  of  the 

game. 

These  three  methods  of  critique,  if  correlated,  will  give  an  excellent 

"replay"  of  the  simulation.  Any  communication  with  the  program  other  than 

listed  above  will  be  accomplished  as  described  by  Leach  and  Perrella  [5]. 

TIMING 

The  timing  of  this  simulation  is  done  by  means  of  a  time  loop  block. 
In  this  block,  the  contents  of  cell  5006B  in  the  CDC  1604- is  tested  and 
stored  as  ICLOCK.  A  variable,  NEXT,  is  generated  as  the  sum  of  ICLOCK 
and  ISTEP,  the  stepping  interval.   ISTEP  is  determined  from  another  variable 
TFACTOR  that  is  equal  to  one  for  real  time.  TF ACTOR  equal  to  three  would, 
for  example,  cause  the  game  to  run  at  three  times  real  time.  When  ICLOCK 
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becomes  greater  than  or  equal  to  NEXT  the  loop  is  exited  and  the  simu- 
lation continues.  This  procedure  is  easily  followed  in  Appendices 
II  and  IV. 


37 


6.   CONCLUSIONS  AND  ACKNOWLEDGEMENTS 

The  considerations  made  in  Che  formulation  and  construction  of 
this  simulation  have  made  some  observations  possible.  Only  the  basic 
modular  structure  of  the  simulation  has  been  considered  in  detail.   Each 
individual  modular  block  has  been  designed  in  as  simple  a  form  as 
possible  while  maintaining  some  degree  of  realism.  The  simulation  is 
sound  in  its  general  organization.  New  program  blocks  may  be  substituted 
making  the  simulation  as  realistic  as  desired.   It  is  hoped  that  this 
simulation  will  be  played  with  more  sophisticated  models  and  on-line 
equipment  of  greater  capacity  so  that  doctrine  and  tactics  in  the  area  of 
self-inflicted  nuclear  environment  may  be  explored. 

It  was  found  that  the  general  purpose  computer  language  was  com- 
pletely satisfactory  for  the  construction  of  this  simulation.  The  mod- 
ularization technique  made  the  logical  organization  of  the  simulation 
straight  forward  and  is  recommended  for  use  in  future  simulations. 

The  author  wishes  to  express  his  appreciation  to  Professor 
Mitchell  L.  Cotton  and  Professor  Alvin  F.  Andrus  for  their  aid  and 
encouragement  in  the  preparation  of  this  thesis.  In  addition,  the  author 
would  like  to  thank  Miss  Patricia  Hoang  for  her  assistance  in  programming 
the  linking  subroutines  and  LCDR  Richard  E.  DeWinter  for  his  comments 
in  proofreading  the  Manuscript. 
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APPENDIX  I 


LIST  OF  VARIABLES 
This  Appendix  contains  a  listing  of  the  variables  used  in  the 
simulation  MULNUCl,  arranged  in  alphabetical  order. 
Variable    Definition 


AROCMAX 


AROCMIN 


BNPTS 

CENTERB 

CLOUDR(I) 
CLOUDX(I) 

CLOUDY (I) 

CONTB(I) 

CONTC(I) 

CONTR(I) 
COOTS  (I) 

COOTX(I) 

COOTY(I) 


Maximum  range  of  the  ASROC,  in  yards.  This  is  an 
input  parameter,  set  equal  to  9,000  yards  by  the  inputs 
block  and  can  be  changed  in  the  change  inputs  block. 

Minimum  safe  range  for  the  ASROC,  in  yards.  This  value 
is  computed  in  the  initializing  block,  and  is  a  function 
of  warhead  size. 

Last  bearing  the  submarine  held  of  the  nearest  destroyer, 
in  degrees  true. 

A  dummy  variable  used  in  the  contact  tracking  block 
for  the  determination  of  contact  speed. 

The  course  such  that  the  submarine  will  split  the 
channel  between  the  destroyers,  in  degrees  true. 

Radius  of  the  i   cloud  of  radiation,  in  yards. 

X-coordinate  of  the  itn  cloud  of  radiation  with  respect 
to  the  main  coordinate  system,  in  yards. 

Y-coordinate  of  the  i   cloud  of  radiation  with  respect 
to  the  main  coordinate  system,  in  yards. 

th 

Bearing  of  sonar  contact  of  i   destroyer,  in  degrees 
true. 

Contact  course,  as  computed  by  the  contact  tracking 
block,  in  degrees  true. 

Sonar  range  to  contact,  as  measured  by  the  i   destroyer. 

Contact  speed  as  computed  by  the  contact  tracking  model, 
in  knots. 

X-coordinate  of  the  i   destroyer's  contact,  in  yards, 
with  respect  co  the  main  coordinate  system. 

Y-coordinate  of  the  i   destroyer  s  contact,  in  yards, 
with  respect  to  the  main  coordinate  system. 
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Variable    Definition 

D  Difference  between  the  true  bearing  to  the  submarine 

from  the  destroyer  and  the  submarine's  true  course, 
in  degrees.  This  is  used  in  the  determination  of 
doppler. 

DA         Absolute  value  of  D. 

DAMAGE      Percent  damage  to  the  submarine  from  the  current 
detonation. 

DAMAGET  Cumulative  damage  to  the  submarine,  in  percent. 

DDC(I)  Course  of  the  i   destroyer,  in  degrees  true. 

DDS(I)  Speed  of  the  i   destroyer,  in  knots. 

DDSMAX(I)  Maximum  speed  available  to  the  i  J  destroyer,  in  knots. 

DDX(I)      X-coordinate  of  the  i   destroyer,  in  yards,  with 
respect  to  the  main  coordinate  system. 

DDY(I)      Y-coordinate  of  the  i   destroyer,  in  yards,  with 
respect  to  the  main  coordinate  system. 

DETR       Detection  range,  in  yards.  This  is  a  random  variable 
with  mean  DETRM  and  normally  distributed  with  sigma 
of  .3  times  DETRM. 

DETRM      Mean  detection  range,  in  yards.  A  function  of  GRAD 
SUBD. 

DISTC(I,J)  Distance  of  the  i   destroyer  from  the  center  of  the 
jtn  cloud  of  radiation. 

DISTP(I,J)  Distance  of  the  i   destroyer  from  the  center  of  the 
jth  pool  of  radiation. 

DOB        Depth  of  burst,  in  feet,  of  ASROC  warhead. 

DR         Advance  of  the  destroyer  considered,  in  yards,  between 
marks  as  considered  by  the  tracking  model. 

DTIME      Total  time  delay  from  the  firing  of  an  ASROC  and  the 
detonation,  in  seconds.  This  includes  time  of  flight 
and  time  of  sinking. 

DUMMY      A  dummy  variable  used  throughout  the  program. 

DX         East-west  advance  of  the  destroyer  considered,  in  yards, 
between  marks  as  considered  by  the  tracking  model. 
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Variable    Definition 

DY         North-south  advance  of  the  destroyer  considered,  in 
yards,  between  marks  as  considered  by  the  tracking 
model. 

ESR        Effective  sonar  range,  in  yards.  Computed  by  the 
initializing  block. 

GRAD       Water  temperature  gradient,  in  degrees  per  hundred 
feet  of  depth. 

GTIME       Game  time,  in  minutes.  At  the  start  of  the  game 
GTIME  is  zero. 

GZX  X-coordinate  of  ground  zero  fo^  the  detonation 
considered,  in  yards,  with  respect  to  the  main 
coordinate  system. 

GZY  Y-coordinate  of  ground  zero  for  the  detonation 
considered,  in  yards,  with  respect  to  the  main 
coordinate  system. 

HULL       Submarine  hull  thickness,  in  inches  of  steel. 

IASROC      Hull  number  of  the  destroyer  that  fired  the  last 

ASROC.  If  no  ASROC  has  been  fired,  IASROC  is  zero. 

ICLASS(I)  Sonar  classification  of  the  itn  destroyer's  contact. 
0  -  no  contact,  1  -  possible  submarine,  2  -  probable 
submarine. 

ICLOCK  Contents  of  location  5006B  in  the  CDC  1604.  The  CDC 
1604  steps  this  cell  once  every  second. 

ICIRCX(I)   A  dummy  vector  of  X-coordinates  of  points  to  transmit 
circles  to  the  display. 

ICIRCY(I)   A  dummy  vector  of  Y-coordinates  of  points  to  transmit 
circles  to  the  display. 

ICONTB(I)  Fixed  point  version  of  CONTB(I) ,  used  to  transmit  to 
the  display. 

ICONTC(I)  Fixed  point  version  of  CONTC(I) ,  used  to  transmit  to 
the  display. 

ICRIT#      A  series  of  critique  indicators.   0  -  pass. 

1  1  -  sonar  contact 

2  1  -  ASROC  fired 
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Variable 


Definition 


ICRIT# 
3 
4 
5 
6 
7 
8 
9 
11 
12 

13 

14 
15 
16 
17 
IDDC(I) 

IDDS(I) 

IDDX(I) 

IDDY(I) 

IDEEP 


-  torpedo  fired 

-  ASROC  misfire 

-  target  in  range 

-  sub  sunk 

-  destroyer  sunk 

-  target  too  close  to  shoot  at 

-  target  out  of  ASROC  range 

-  game  is  a  draw,  submarine  escaped 

-  submarine  wins,  destroyer  1  sunk  with  submarine 
escaping 

-  submarine  wins,  destroyer  2  sunk  with  submarine 
escaping 

-  submarine  wins,  both  destroyers  sunk 

-  destroyers  win,  submarine  sunk  by  destroyer  1 

-  destroyers  win,  submarine  sunk  by  destroyer  2 

-  submarine  wins  by  transiting  between  the  destroyers 

ixed  point  version  of  DDC(I),  used  to  transmit  to 

he  display. 

Fixed  point  version  of  DDS(I) ,   used  to  transmit  to  the 
display. 

Fixed  point  version  of  DDX(I)  ,   used  to  transmit  to  the 
display. 

Fixed  point  version  of  DDY(I),   used  to  transmit  to  the 
display. 

An  index  used  to  indicate  the  classification  of  nuclear 
burst.    1  -  very  shallow,   2  -  shallow,   3  -  deep,   4  - 
very  deep. 

IDOPLER(I)     An  index  used  to  indicate  the  doppler  of  the  i 

destroyer's  contact.     0  -  no  doppler,    1  -  up  doppler, 
2  -  down  doppler. 
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Variable    Definition 

IEND       An  index  used  to  indicate  game  over.  0  -  game  not  over, 

1  -  game  over. 

IFIRE      An  index  used  to  indicate  if  the  submarine  has  fired  a 
torpedo.   0  -  torpedo  not  active,  1  -  torpedo  still 
active. 

ILOGIC      A  logic  index  used  in  the  submarine  model. 

IONE       The  lowest  destroyer  hull  number.  Equal  to  1  at  the 

beginning  of  the  game.   If  destroyer  number  one  is  sunk 
the  IONE  is  equal  to  two. 

INITIAL     An  index  of  the  initial  situation.   0  -  submarine  at 
the  origin,  1  -  submarine  randomly  distributed  in 
the  upper  half  of  the  playing  area. 

IR         Radius  of  display  area,  in  thousands  of  yards. 

IRANDOM    A  random  number  selected  to  initialize  the  random 

number  generator.  This  number  must  be  an  odd  integer 
in  the  interval  1  to  67108863. 

ISHOOTl     A  flag  used  to  activate  the  firing  sequence  for 
destroyer  1.  0  -  do  not  shoot,  1  -  shoot. 

ISH00T2     A  flag  used  to  activate  the  firing  sequence  for 
destroyer  2.  0  -  do  not  shoot,  2  -  shoot. 

ISOL(I)     A  progressive  index  of  the  quality  of  the  firing 

solution  the  i*-*1  destroyer  has  on  its  target.  0  -  no 
solution,  thru  5  -  best  solution. 

ISS        Sea  state. 

ISTEP      Fixed  point  version  of  TSTEP. 

ISTRAT      Submarine  basic  strategy  when  the  submarine  is  controlled 
by  the  computer.  0  -  run  for  it,  1  -  go  up  the  middle, 

2  -  end  run. 

ISUBC      Fixed  point  version  of  SUBC. 

ITURN      Logical  index  that  records  the  submarine's  initial 
turn.  0  -  left  turn,  1  -  right  turn. 

1X0        X-coordinate  of  center  of  displayed  area,  in  hundreds 
of  yards  with  respect  to  the  main  coordinate  system. 

IXDD(I,J)   A  dynamic  table  of  track  data,  recording  the  X-coordinate 
of  the  jtn  mark  of  the  itn  destroyer.  The  maximum  value 
of  j  is  eight.   Entries  are  in  yards,  with  respect  to  the 
main  coordinate  system. 
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Variable 


Definition 


IXDDl(I) 
IXDD2(I) 
IXSUB(I,J) 


IXSUBl(I) 
IXSUB2(I) 
IYO 

IYDD(I,J) 

IYDDl(I) 
IYDD2(I) 
IYSUB(I,J) 

IYSUBI(I) 
IYSUB2(I) 
IWINDD 

MARKS  (I) 

N 
NCONBER 

NEXT 
NMARKS(I) 


This  variable  is  equal  to  IXDD(1,  I) . 

This  variable  is  equal  to  IXDD(2,  I). 

A  dynamic  table  of  track  data,  recording  the  X- 
coordinate  of  the  jth  mark  of  the  itn  destroyer's 
sonar  contact.  The  maximum  value  of  j  is  eight. 
Entries  are  in  yards,  with  respect  to  the  main  coordinate 
system. 

This  variable  is  equal  to  IXSUB(1,I). 

This  variable  is  equal  to  IXSUB(2,I). 

Y-coordinate  of  center  of  displayed  area,  in  hundreds 
of  yards  with  respect  to  the  main  coordinate  system. 

A  dynamic  table  of  track  data,  recording  the  Y-coordinate 
of  the  jfch  mark  of  the  ith  destroyer.  The  maximum 
value  of  j  is  eight.   Entries  are  in  yards,  with  respect 
to  the  main  coordinate  system. 

This  variable  is  equal  to  IYDD(]»I)  . 

This  variable  is  equal  to  IYDD(2,I). 

A  dynaoic  table  of  track  data,  recording  the  Y-coordinate 
of  the  jtn  mark  of  the  itn  destroyer's  sonar  contact. 
The  maximum  value  of  j  is  eight.  Entries  are  in  yards, 
with  respect  to  the  main  coordinate  system. 

This  variable  is  equal  to  IYSUB(1,  I). 

This  variable  is  equal  to  IYSUB(2,  I). 

Fixed  point  version  of  WINDD,  used  to  transmit  data  to 
the  display. 

Number  of  continuous  marks,  up  to  five,  the  i 
destroyer  has  on  his  sonar  contact. 

A  dummy  variable. 

Number  of  constant  bearing  the  submarine  has  on  the 
nearest  destroyer. 

Time  of  the  next  cycle,  in  seconds  of  computer  time. 

Same  as  MARKS(I)  except  NMARKS(I)  does  not  stop  at  five. 
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Variable 


Definition 


NOSHOOT 

NPTS 

NPTSl 
NPTS  2 
NRDD 

NSHOTS 
PHIT 

POOLR(I) 
POOLX(I) 

POOLY(I) 

R 

REL 

RADDOSE(I) 

RADRATE(I) 

RANDOM 

RANGE 

RLETHAL 
SAFETY 

SB(I) 


Index  to  limit  destroyers  to  one  active  weapon  at  a 
time.   0  -  there  are  no  active  weapons,  alright  to 
shoot,  1  -  there  is  an  active  weapon,  cannot  shoot. 

Index  used  to  control  display  data  until  there  are 
eight  point  available  on  destroyer  tracks. 

A  dummy  variable. 

A  dummy  variable. 

The  highest  destroyer  hull  number.  Equal  to  2  at  the 
beginning  of  the  game.   If  destroyer  number  two  is 
sunk  then  NRDD  is  equal  to  one. 

Total  number  of  ASROC' s  fired  during  the  game. 

A  random  variable  used  to  determine  if  a  torpedo,  fired 
by  the  submarine,  hit  the  destroyer. 

Radius  of  the  i   pool  of  radiation,  in  yards. 

X-coordinate  of  the  i"1  pool  of  radiation,  in  yards, 
with  respect  to  the  main  coordinate  system. 

Y-coordinate  of  the  i   pool  of  radiation,  in  yards, 
with  respect  to  the  main  coordinate  system. 

Radius  of  the  displayed  area,  in  yards. 

A  random  variable  used  to  determine  ASROC  reliability. 

Total  radiation  dose  the  i   destroyer  has  been  exposed 
to  during  the  game,  in  roentgens. 

The  rate  at  which  the  i   destroyer  is  receiving 
radiation,  in  roentgens/hr. 

The  random  variable  used  to  link  the  various  random 
generators. 

A  dummy  variable  used  as  a  temporary  storage  in  range 
calculations. 

The  lethal  range  of  the  ASROC  warhead,  in  yards. 

Safety  factor  used  to  compute  submarine  maximum 
operating  depth  as  a  function  of  the  crush  depth  of 
the  hull. 

Screw  beat  count  of  the  l"1  destroyer  as  measured  by  the 
submarine. 
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Variable    Definition 

SIGMA      A  dummy  variable  used  in  the  calling  of  normally 
distributed  random  variables. 

SR         Sink  rate  of  the  ASROC  warhead,  in  feet  per  second. 

SSB(I)      The  actual  beating  of  the  submarine  from  the  i 
destroyer,  in  degrees  true. 

STRESS      Yield  strength  of  the  steel  used  in  the  submarine 

hull,  measured  in  thousands  of  pounds  per  square  inch. 

SUBC       The  true  submarine  course,  in  degrees  true. 

SUBD        The  actual  submarine  depth  in  feet. 

SUBDMAX     The  maximum  allowable  operational  depth  of  the  submarine, 
in  feet. 

SUBS        The  actual  speed  of  the  submarine,  in  knots. 

SUBSMAX     The  maximum  speed  available  to  the  submarine,  in  knots. 

SUBX       The  actual  submarine  X-coordinate  with  respect  to  the 
main  coordinate  system,  in  yards. 

SUBY       The  actual  submarine  Y-coordinate  with  respect  to  the 
main  coordinate  system,  in  yards. 

SUM        A  dummy  variable  used  in  the  computation  of  radiation 
dose. 

TBURST(I)   Time  of  detonation  of  the  ith  ASROC  warhead,  measured 
in  game  time. 

TEMP       A  dummy  variable. 

TEMPI       A  dummy  variable. 

TEMP2       A  dummy  variable. 

TEMP3      A  dummy  variable. 

TEMP4       A  dummy  variable. 

A  dummy  variable. 


TEMP5 

TEMPCR(I,J)  The  radiation  rate  received  by  the  iu"  destroyer  from 
ith 


th 


the  jzn   cloud  of  radiation,  in  roentgens  per  hour. 

TEMPPR(I,J)  The  radiation  rate  received  by  the  itn  destroyer  from 
the  jtn  cloud  of  radiation,  in  roentgens  per  hour. 
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Variable    Definition 

TFACTOR     A  time  factor,  1  -  real  time,  2  -  double  time,  .5  - 
half  time. 

THERMO      Depth  of  the  thermocline,  in  feet. 

THETA      A  dummy  variable  used  in  the  computation  of  courses  and 
bearings. 

TINTER      Time  to  intercept  of  a  torpedo  fired  by  the  submarine, 
in  minutes. 

TLOGIC      A  time  storage  point  used  to  control  timed  logic. 

TOB        Time  of  burst,  this  includes  time  of  flight  and  sink 

time  of  the  ASROC  warhead,  measured  in  minutes  of  game 
time. 

TOF        Time  of  flight  of  the  ASROC  warhead,  in  seconds. 

TOS        Sink  time  of  the  ASROC  warhead,  in  seconds. 

TSTEP      The  time  step  for  each  cycle  of  the  game,  in  seconds. 

TXDD(I,J)  A  temporary  storage  of  the  j  past  position  X-coordinate 
of  the  ifch  destroyer  used  in  preparing  track  data  for  the 
display,  measured  in  yards. 

TYDD(I,J)  A  temporary  storage  of  the  j*-*1  past  position  Y-coordinate 
of  the  itn  destroyer  used  in  preparing  track  data  for  the 
display,  measured  in  yards. 

TXSUB(I,J)  A  temporary  storage  of  the  j   past  position  X-coordinate 
of  the  itn  destroyer's  sonar  contact  used  in  preparing 
track  data  for  the  display,  measured  in  yards. 

TYSUB(I,J)  A  temporary  storage  of  the  j"1  past  position  Y-coordinate 
of  the  itn  destroyer's  sonar  contact  used  in  preparing 
track  data  for  the  display,  measured  in  yards. 


VEL 

Wl 

W2 

WINDD 
WINDV 


The  relative  velocity  of  the  torpedo  fired  by  the  sub- 
marine, with  respect  to  the  destroyer  target,  in  knots. 

A  dummy  variable  used  to  determine  the  destroyer  that 
is  closest  to  the  submarine. 

A  dummy  variable  used  to  determine  the  destroyer  that 
is  closest  to  the  submarine. 

Wind  direction  in  degrees  true. 

Wind  velocity,  in  knots. 
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Variable 
X(I,J) 

XDD(I,J) 

xo 

XTEMP 
Y(I,J) 


Definition 


The  X-coordinate  of  the  i   destroyer's  j   continuous 
mark  on  its  sonar  contact,  in  yards,  with  respect  to 
the  main  coordinate  system. 

The  X-coordinate  of  the  itn  destroyer  j-1  steps  ago, 
in  yards,  with  respect  to  the  main  coordinate  system. 

The  X-coordinate  of  the  center  of  the  displayed  area, 
in  yards,  with  respect  to  the  main  coordinate  system. 

A  dummy  variable  used  for  temporary  storage  of  X-coordinates. 

The  Y-coordinate  of  the  ith  destroyer's  jth  continuous 
mark  on  its  sonar  contact,  in  yards,  with  respect  to  the 
main  coordinate  system. 

th 

YDD(I,J)    The  Y-coordinate  of  the  1   destroyer  j-1  steps  ago, 
in  yards,  with  respect  to  the  main  coordinate  system. 

YIELD      The  yield  of  the  ASROC  warhead  in  kilotons. 

YO         The  Y-coordinate  of  the  center  of  the  displayed  area, 
in  yards,  with  respect  to  the  main  coordinate  system. 

ZZX        A  vector  of  destroyer  and  submarine  X-coordinates  used 
by  the  graph  plot  routine  in  Critique  I. 

ZZY        A  vector  of  destroyer  and  submarine  Y-coordinates  used 
by  the  graph  plot  routine  in  Critique  I. 
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APPENDIX  II 
LOGIC  FLOW  DIAGRAMS 

This  Appendix  contains  a  series  of  logic  flow  diagrams  as  listed 
below. 

Executive  Control 
Plot  Generator 
Interactions 
Radiation  Model 
Submarine  Logic  Model 
Run  for  It 
Up  the  Middle 
End  Run 

Submarine  Team  Control 
Sonar  Contact  Model 
Contact  Tracking  Model 
Weapon  Firing  Model 
Evaluation  Model 
Critique  I 
Critique  II 
It  will  be  noted  that  the  statement  number  from  the  program  listing  is 
shown  in  the  upper  right  hand  corner  of  each  symbol  in  the  block  diagram. 
This  should  aid  in  correlating  this  Appendix,  with  Appendix  IV. 
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FLOW  CHART   SYMBOLS 


CD 


A  connector  or  terminal. 


An  offpage  connector. 


A  predefined  process  or  module/ subroutine. 

A  more  detailed  flow  chart  of  this  subroutine  is  also 

included. 


Input/output  other  than  display. 


Decision. 


Processing,  annotation. 


Display. 
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APPENDIX  III 
SUBROUTINES  AND  CDC  160  EXECUTIVE  ROUTINE 

This  Appendix  contains  as  explanation  of  the  subroutines  used  in 
the  main  program  and  the  CDC  160  executive  routine  used  to  connect  the 
CDC  1604  with  the  DD  65  display.  A  listing  of  these  subroutines  can  be 
found  in  Appendix  IV. 

SUBROUTINE  RANVAR 
Subroutine  RANVAR  is  used  as  the  basic  random  number  generator.  It 
generates  random  floating  point  numbers  in  the  interval  zero  to  one  that 
are  distributed  uniformly  in  that  interval.  The  random  variables  used 
in  the  main  program  are  called  from  either  UNIFORM,  NORMAL  or  ERROR 
subroutines  which  in  turn  call  RANVAR  for  input.  The  generator  is  a 
simple  fixed  point  division  utilizing  the  remainder  from  the  Q  register 
as  the  random  number.  This  number  is  then  mapped  into  the  zero  to  one 
interval.  Only  one  input  is  required  to  initialize  this  generator, 
namely  IRANDOM. 

SUBROUTINE  UNIFORM 
Subroutine  UNIFORM  is  a  three  argument  subroutine  used  to  generate 
uniformly  distributed  random  numbers  in  any  interval.  The  arguments  are: 

1.  The  center  of  the  interval. 

2.  The  half  width  of  the  interval. 

3.  The  output  random  number. 

An  example  would  be  UNIFORM  (5.0,  2.0,  SUBS).  This  call  would  yield  the 
variable  SUBS,  submarine  speed,  uniformly  distributed  in  the  interval 
with  center  at  five  knots  and  plus  or  minus  two  knots. 
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SUBROUTINE  NORMAL 
Subroutine  NORMAL  is  a  three  argument  subroutine  used  to  generate 
normally  distributed  random  numbers  with  any  mean  and  standard  deviation. 
This  normal  distribution  is  generated  by  means  of  the  sum  of  identically 
distributed  (uniform)  random  variables.  Twelve  uniform  numbers  are  used 
because: 

1.  The  truncation  is  not  significant. 

2.  Twelve  reduces  the  formula  to  a  summation  and  no  division  is 
required. 

The  arguments  to  this  subroutine  are: 

1.  The  mean  of  the  distribution. 

2.  The  signia  of  the  distribution. 

3.  The  output  random  number. 

SUBROUTINE  ERROR 
Subroutine  ERROR  is  a  five  argument  subroutine  used  to  generate 
circular  normal  distributed  random  variables.  This  subroutine  is  used 
to  determine  the  true  fall  of  shot  given  the  aiming  point  and  CEP  or 
sigma  of  the  distribution  of  fall  or  shot.  The  arguments  are: 

1.  The  X-coordinate  of  the  aim  point. 

2.  The  Y-coordinate  of  the  aim  point. 

3.  The  sigma  of  the  circular  normal  distribution. 

4.  The  X-coordinate  of  the  fall  of  shot. 

5.  The  Y-coordinate  of  the  fall  of  shot. 

SUBROUTINE  CIRCLE 
Subroutine  CIRCLE  is  used  to  generate  a  series  of  points,  24  in 
number,  every  15  degrees  around  the  perimeter  of  a  circle  of  predetermined 
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center  and  radius.  The  circle  is  then  used  in  the  formation  of  the 
circular  clouds  and  pools  of  radiation  displayed  to  the  destroyer  team. 
The  arguments  are: 

1.  The  X-coordinate  of  the  center  of  the  circle  with  respect  to 
the  main  coordinate  syste. 

2.  The  Y-coordinate  of  the  center  of  the  circle  with  respect  to 
the  main  coordinate  system. 

3.  The  radius  of  the  circle  with  respect  to  the  main  coordinate 
system. 

4.  The  X-coordinate  of  the  center  of  the  display  with  respect  to 
the  main  coordinate  system. 

5.  The  Y-coordinate  of  the  center  of  the  display  with  respect  to 
the  main  coordinate  system. 

6.  The  radius  of  the  displayed  area. 

7.  A  vector  of  X-coordinates  of  the  24  points  in  the  circle  with 
respect  to  the  display  coordinate  system. 

8.  A  vector  of  Y-coordinates  of  the  24  points  in  the  circle  with 
respect  to  the  display  coordinate  system. 

SUBROUTINE  DCIRCLE 
DCIRCLE  is  a  subroutine  used  to  transmit  the  circle  coordinates, 
generated  in  subroutine  CIRCLE,  to  the  CDC  160  from  the  CDC  1604.  The 
arguments  to  this  subroutine  are: 

1.  ITRKNO  -  the  track  of  circle  number  by  which  this  particular 
circle  can  be  designated. 

2.  CHAR  -  a  single  letter  or  number  in  hollerith  form  that  is  to 
be  displayed  as  one  of  the  points  in  the  circle 

3.  In  the  program  MULNUCl  the  letter  c  and  p  are  used  to  dis- 
tinguish between  clouds  and  pools  of  radiation. 
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4.  NUMPTS  -  the  number  of  points  in  the  circle. 

5.  IX  -  a  vector  of  X-coordinates  for  the  circle. 

6.  IY  -  a  vector  of  Y-coordinates  for  the  circle. 

SUBROUTINE  DTRACK 
Subroutine  DTRACK  is  used  to  transmit  track  data  to  the  CDC  160 
from  the  CDC  1604.  This  subroutine  will  send  tracks  of  up  to  eight  points 
to  the  CDC  160.  The  arguments  to  this  subroutine  are  the  same  as  those 
in  DCIRCLE,  in  fact  the  two  routines  are  identical  with  the  exception  of 
the  allowable  number  of  points. 

SUBROUTINE  DSTATUS 
DSTATUS  is  a  subroutine  used  to  transmit  information  to  the  windows 
described  in  the  section  on  Man  Machine  Interface.  This  information  can 
be  in  two  forms: 

1.  Program  variables  in  fixed  point  form. 

2.  Messages  in  hollerith  form. 
The  arguments  are: 

1.  ITYPE  -  zero  represents  a  numerical  program  variable  is  to  be 
sent,  while  one  indicates  that  an  eight  hollerith  character  word  is  to  be 
sent. 

2.  NWIND  -  the  number  of  the  window  to  which  the  variable  or  word 
is  to  be  sent. 

3.  IW  -  the  field  width  of  the  variable. 

4.  INAME  -  the  name  of  the  variable  to  be  transmitted. 

5.  IX  -  the  X-coordinate  of  the  lower  left  hand  corner  of  the 
window. 

6.  IY  -  the  Y-coordinate  of  the  lower  left  hand  corner  of  the 

window. 
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It  should  be  noted  that  windows  are  128  display  coordinates  units  long. 
If  a  message  in  the  form  of  words  is  to  be  sent  to  the  display  and  it  is 
longer  than  eight  letters  (the  length  of  one  window)  it  can  be  sent  by 
means  of  more  than  one  window.  These  windows  should  be  displaced  by  128 
units  in  the  X  direction,  thus  the  windows  may  form  a  continuous  word  of 
more  than  eight  characters. 

SUBROUTINE  PARAMS 
PARAMS  is  an  eight  argument  subroutine  used  to  query  the  CDC  160 
as  to  the  contents  of  eight  selected  windows.  This  subroutine  allows 
the  main  program,  in  the  CDC  1604,  to  enter  changes  that  have  been  made 
to  the  windows  of  the  display  by  the  player.  This  is  the  only  method  the 
player  has  of  communicating  with  the  program  without  interrupting  the 
play.  The  eight  arguments  to  this  subroutine  are: 

lA  -  the  contents  of  window  5 

2A  -  the  contents  of  window  6 

3A  -  the  contents  of  window  7 

4A  -  the  contents  of  window  8 

5A  -  the  contents  of  window  9 

6A  -  the  contents  of  window  9 

7A  -  the  contents  of  window  16 

8A  -  the  Burner ic  value  of  a  location  in  the  CDC  160  that  is 
controlled  by  the  DDl  FIRE  and  DD2  FIRE  buttons. 
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APPENDIX  IV 
PROGRAM  LISTING 
This  Appendix  contains  the  computer  program  listing  of  the  simu- 
lation MULNUCl.  The  program  blocks  are  in  numerical  order  and  the  logic 
can  be  followed  by  cross  referencing  this  Appendix  with  Appendix  I.   In 
this  program  all  classified  input  parameters  have  been  assigned  fictitious 
values  so  that  the  program,  as  presented  in  this  thesis,  could  remain 
unclassified. 
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